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In the last years there has been a renewed interest for zinc oxide semiconductor, mainly triggered by its 
prospects in optoelectronic applications. In particular, zinc oxide thin films are being widely used for 
photovoltaic applications, in which the determination of the electrical conductivity is of great importance. 
Being an intrinsically doped material, the quantification of its doping concentration has always been 
challenging. Here we show how to probe the charge carrier density of zinc oxide thin films by Scanning 
Kelvin Probe Microscopy, a technique that allows measuring the contact potential difference between the tip 
and the sample surface with high spatial resolution. A simple electronic energy model is used for correlating 
the contact potential difference with the doping concentration in the material. Limitations of this technique 
are discussed in details and some experimental solutions are proposed. Two-dimensional doping 
concentration images acquired on radio frequency-sputtered intrinsic zinc oxide thin films with different 
thickness and deposited under different conditions are reported. We show that results inferred with this 
technique are in accordance with carrier concentration expected for zinc oxide thin films deposited under 
different conditions and obtained from resistivity and mobility measurements. 



In the last century, Zinc Oxide (ZnO) attracted interest within the scientific community. Defined as the future 
material, ZnO has been widely studied since 1935. Lately a renewed interest for this material had risen. This is 
mainly due to reports ofp-type conduction\ ferromagnetic behavior^ and improvements in epitaxial growth 
technologies that make this material suitable for the fabrication of optoelectronic^, piezoelectric^ and spintronic 
devices^. Extensive reviews on ZnO applications are available in the literature^'^. Nano- crystalline ZnO thin films 
are being widely studied, especially for photovoltaic applications^'^. Due to its inherent transparency in the visible 
range, ZnO is used as transparent conductive oxide (TCO) on top of thin film solar celP°'^\ As an example, in 
inverted bulk-heteroj unction organic solar cells^^'^^ n-doped ZnO can be used as the electron collector/hole 
blocking buffer layer^^. 

ZnO thin films can be either deposited^^ or solution processed^^. Determination of the doping concentration of 
ZnO films is a fundamental step for finely tuning ZnO deposition parameters (e.g. for epitaxial growth) and ZnO 
nanoparticles synthesis (for solution processing). A method able to determine thin films ZnO doping profiles at 
the nanoscale would allow optimizing its performance in several applications. 

Although different techniques to measure the doping concentration in semiconductors have been 
developed^^"^^, the determination of this property in ZnO thin films still remains challenging^^. This is due to 
the fact that ZnO is intrinsically n-doped, presumably by Zn interstitials, O vacancies and H, that affect the 
position of the Fermi level within the bandgap^". In this view, ordinary method such as Secondary Ion-Mass 
Spectroscopy (SIMS) and Electron Holography (EH) can just give an insight of the possible cause of the self- 
doping effect^ without giving accurate quantitative results in terms of charge carrier concentration. 
Furthermore, these two techniques are destructive and therefore their versatility is limited. Four Probe 
Method is the most widespread technique to measure the resistivity of thin films^^. From this measurement, 
given the carriers mobility, the effective doping concentration of the material can be recovered. However, several 
experimental issues have to be taken into account when using this technique. The contact between a metal and a 
semiconductor can form a Schottky diode^^ rather than an ohmic contact. To obtain reliable results, adjustments 
of the drive current are needed depending of the resistivity of the sample. Due to the reduced thickness of thin 
films, the measurement can be affected by the presence of the substrate, leading to erroneous results. In this view, a 



SCIENTIFIC REPORT! | 4 : 4203 | DOI: 1 0.1 038/srep04203 



1 



non-contact technique able to measure the position of the Fermi level 
and to relate it with the effective doping concentration would rep- 
resent an adequate solution for the aforementioned issues. 

Scanning Kelvin Probe Microscopy (SKPM), a modified version of 
Atomic Force Microscopy (AFM), is a non-destructive non-contact 
surface technique that allows imaging two-dimensional profiles of 
contact potential difference (Vcpd) i-e. the difference in the work 
functions of the tip and the sample^^"^^. This is true in absence of 
trapped charges^^ or adsorbates on the surface. It is well known that 
in semiconductors the work function depends on the dopant types 
and concentrations^^. This is because they affect the Fermi level 
position within the bandgap. 

Different attempts to correlate SKPM measurements with doping 
concentration in semiconductors have been reported in literature 
since the technique has been proposed^^'^°. Henning^^ Shin^^ and 
Hochwitz^^ reported a doping concentration two-dimensional image 
in silicon microstructures and integrated circuits. Due to its inherent 
nanometric spatial resolution^^, SKPM has been successively used 
to probe non-uniform doping distribution at the nanoscale. 
Semenikhin^^ reported evidence of local doping inhomogeneity on 
conducting polybithiophene, while Koren^^ measured the doping 
distribution along single phosphorous-doped silicon nanowires. 
Although compared with Scanning Capacitance Microscopy (SCP) 
this technique presents some limitations in terms of accuracy^^'^^, 
previously mentioned and other reports^^ demonstrate that SKPM 
can be used to quantitatively measure the doping concentration in 
semiconductors with high spatial resolution^^'^^'^^. Attempts to cor- 
relate SKPM measurements with charge carriers concentration in 
ZnO have been reported^°'^\ However, a clear and complete form- 
alism is still lacking. 

Measuring the doping concentration of a sample without touching 
its surface represents a considerable advantage, especially for thin 
films. This avoids damaging the surface and does not introduce extra 
effects due to Schottky barrier formation. Compared to Four Probe 
Method, SKPM can be used to measure the doping profile of semi- 
conducting thin films deposited on metal^\ This would prevent from 
the need to deposit a sample on an insulating substrate just to mea- 
sure its resistivity. Being a surface technique, SKPM is able to mea- 
sure carrier concentration on top of the sample. Here, the presence of 
charges can affect the measurement. To minimize this effect, the 
sample is scanned in tapping mode while both the tip and the back 
contact are connected to the ground, before performing the actual 
experiments. By creating a closed circuit, the surface charges flows 
towards the electric ground, causing the surface to discharge. This 
technique has been successfully applied by Jespersen et al.^^. While 
doping concentration is normally associated with bulk intrinsic or 
extrinsic doping, due to their small thickness, surface defects in thin 



films have a significant contribution on the overall sample doping^^. 
This has been experimentally demonstrated by photovoltage mea- 
surements performed on ZnO nano-"mounds"^^. Since SKPM is able 
to measure the presence of surface defects^^ it represents a suitable 
technique for imaging two-dimensional doping concentration in 
ZnO thin films. 

In this paper, two-dimensional charge carrier concentration 
images acquired on RF- sputtered intrinsic ZnO thin films are 
reported. A simple energy band model that correlates work function 
difference measurements with the effective doping concentration is 
presented. ZnO samples with different thickness (from 10 to 
300 nm) are investigated. To further validate the method, doping 
concentrations of ZnO samples deposited under different conditions 
are compared. We show that results inferred with this technique are 
in accordance with carrier concentration obtained with Four Probe 
and Hall measurements. 

SKPM is an AFM based technique developed by Nonnemacher 
et al. in 1991^^'^°. As other electrostatic force microscopy tech- 
niques^^'^^, it measures the contact potential difference Vcpd between 
a conductive tip and a sample. The Vcpd is defined as: 



Vr 



sample 



(1) 



where q is the electronic charge, while (l)tip and (j) sample are the work 
functions of the tip and the sample respectively. Compared to 
Ultraviolet Photoelectron Spectroscopy (UPS), which measures the 
kinetic energy of spectra of photoelectrons emitted under ultraviolet 
excitation, both techniques can be used to quantify the work function 
of a given material and therefore could be employed to determine the 
charge carrier concentration. However, UPS has a lower spatial reso- 
lution compared to SKPM, being limited by the dimension of the 
ultraviolet light spot used as excitation^^. Thus, in order to quantify 
the charge carrier concentration with high spatial resolution, we rely 
on SKPM. In SKPM, when the tip is brought close to the sample and 
they are both grounded, an electrostatic interaction due to the differ- 
ences in their Fermi energy levels can be probed. Fig. 1 shows the 
energy level diagram of the tip and the sample for different (l)tip and 
^sample- When tip and sample are not electrically connected, their 
vacuum levels are aligned but Fermi energy levels are different 
(Fig. la). Upon electrical contact (for example both grounded), the 
system reaches equilibrium, and Fermi levels line-up through elec- 
tron current flow (Fig. lb). This current flow causes the vacuum 
energy levels to be no longer the same: the consequence of this is 
the formation of a contact potential difference Vcpd that affects the 
interaction between the tip and the sample, causing the formation of 
an electrostatic force between them. As shown in Fig. Ic, this force 




Figure 1 | Electronic energy levels of the tip-sample system for three different cases, (a) tip and sample are separated by a distance d and are not 
electrically connected; (b) tip and sample are electrically connected. In this case the Fermi energy levels are lined-up; (c) an external bias equal to the 
contact potential difference Vcpd is applied to the tip. This causes the tip-sample electrical force to be nullified. 
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can be nullified. This is done indeed applying an external bias V^c = 
VcPD between the tip and the sample. Thus, in this case, the amount 
of applied external bias is exactly equal to the work function differ- 
ence between the two. 

In SKPM the contact potential difference Vcpd between a con- 
ductive tip and a sample is measured by applying to the tip an ac 
voltage (Vac — ^Acsin(ft)O) and a dc voltage (Vdc)- The V^c causes 
the cantilever to oscillate because of the electrical forces between the 
AFM tip and the sample surface, while the Vqc nullifies the oscil- 
lating electrical forces that originated from the Vcpd- The frequency 
CO is chosen equal or close to the resonance frequency of the 
cantilever in order to increase the signal-to-noise ratio. The electro- 
static force component measured at frequency co has the following 
expression^^'^"'^^: 

dc 

Po^=--^ (yDC-VcPD)VAc sm{cot) (2) 

where C is the tip -sample capacitance and z is the tip -sample dis- 
tance. From (2), it is evident that the Vcpd can be measured by 
applying a dc voltage Vdc such that the oscillating amplitude at co 
is nullified. Equation (2) is strictly valid for metallic samples. 
However, it has been demonstrated that a similar expression can 
be derived for semiconductors, provided that the tip-sample capa- 
citance is written as the sum of the air gap and the space- charge 
capacitor^^. 

The Fermi energy level in semiconductors depends on the doping 
concentration^^. Fig. 2 depicts the electronic energy level for a 
tip -semiconductor system when the tip and the back surface of 
the semiconductor are electrically connected. In the case of a non- 
degenerate n-doped semiconductor, the Boltzmann approximation 
for the concentration of electrons (n) reads^^: 




Figure 2 | Electronic energy band line-up for a metallic tip- 
semiconductor system when the tip and the back surface of the 
semiconductor are electrically connected. The contact potential difference 
Vcpd depends on the tip work function on the electron affinity of 
the semiconductor Xy as well as on the difference between the bottom 
energy level of the conduction band Ec and the Fermi level Ep in the 
semiconductor. In this case, the effect of the surface band bending is 
included in the bulk. 



where Nc is the effective density of states in the conduction band, 
and Ep are the conduction band energy level and the Fermi energy 
level, Kb is Boltzmann constant and T is the temperature. If the 
effective donors concentration Cd is significantly higher than the 
intrinsic carriers concentrations W/, the concentration of electrons 
n can be used to estimate Cd by using the following approximation^^: 

Co'^N^-N^'-n (4) 

where the effective donors concentration Cd is defined as the differ- 
ence between the ionized concentration of donors and the 
ionized concentration of acceptors . When the n-doped semi- 
conductor surface is electrically connected to the tip through a back 
contact, the Vcpd can be written as^^: 

V lip-X+{Ep-Ec) 

ycpD= (5; 

where x is the electron affinity of the semiconductor. Analogous 
expressions can be derived for a p-doped semiconducting surface. 
In equation (5) the band bending term is not explicitly included 
and its effect on the electrical conductivity of the surface is taken into 
account by including it in the bulk. Due to their limited thickness, 
indeed, we expect the presence of defect states at the surface to affect 
the charge carrier concentration through the entire film. This 
assumption is supported by Ref. 44, where photovoltage measure- 
ments confirm the presence of electrically- active deep levels due to 
zinc vacancies that affects the surface potential on a scale of hundreds 
of nanometers. 

Finally, by combining equations (3-5), Vcpd images probed by 
SKPM can be mapped into the effective donors concentration Cd, 
which actually corresponds to the effective doping profile: 

C.^iV.exp^ j 

While performing SKPM scans, two main sources of errors must be 
taken into account: stray capacitances from tip geometry, and envir- 
onmentally mediated water layers presence^^. In the ideal case, it is 
generally assumed that the tip-sample capacitance is essentially due 
to the effect of the tip apex^^. However, the contribution from the tip 
cone and the cantilever can be significant^^'^°. Therefore, it is possible 
that measured SKPM values differ from theoretical ones. Gil et al.^^ 
reported that the standard interpretation of data may be incorrect 
because the total force is produced not only by the end of the tip, but 
also by the contributions of the cantilever and the tip cone, which 
senses a different surface potential than the tip. Jacobs et al.^^ demon- 
strated that reasonable Vcpd values can only be obtained by using 
tips with a long cone and by reducing the lift height during the SKPM 
interleave scan. In the present work, tips with a long cone (tip height 
14 ± 4 jim) and a small tip angle (15 ± 5 degrees) are employed. 
Moreover, the surface is scanned at a small tip-sample distance 
(20 nm). This allows minimizing the effect of the tip cone and the 
lever, which would dominate instead at larger tip-sample separa- 
tions. We have quantified indeed that at 20 nm away from the sur- 
face, for our cantilever (ASYELEC 0.1) the tip contributes to more 
than 50% of the total electrostatic interaction and its contribution is 
therefore the most relevant among all. The contributions of the cone 
and the lever can be further reduced by measuring the gradient of the 
force instead of the force itself^. Furthermore, scanning at small tip- 
sample separations allows maintaining a high spatial resolution 
because the smaller is the distance, the smaller is the area of inter- 
action of the AFM probe^^. A small oscillating amplitude (5 nm) is 
used in this case to maintain the system working in the harmonic 
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Figure 3 | (a) Atomic force microscopy topography image of a ZnO thin film deposited by RF sputtering. Small grains with dimensions of tens of 
nanometers can be observed (b) Two-dimensional contact potential difference Vcpd image of the same surface imaged in Fig. 3a. The image is acquired at 
an equilibrium distance of 20 nm with an oscillation of 5 nm. 



regime. This has been verified by looking at the shape of the electrical 
tuning where no distortions in the Lorentzian function (amplitude 
vs. frequency) were observed. When a high value of Vac is employed, 
the tuning curve may present distortions, causing bi-stability in the 
tip-sample interaction^^ Another source of inaccuracy in the mea- 
surement of Vcpd is the presence of water layers or molecular adsor- 
bates on the sample surface. It has been shown that the presence of a 
thin layer of water or water vapor in the environment can cause 
deviations in the measurement of the resistivity of metal-oxide 
films^^. Water molecules are generally chemisorbed and physisorbed 
by the surface^^, and this process increases or decreases the surface 
electron conductivity depending on the type of the semiconductor 
(respectively n ox p type)^^. In the present work, we minimize the 
effect of water layers and adsorbates by preparing the samples by 
cycles of Ar^ bombardment at 250°C, followed by annealing at 
120°C in a vacuum chamber. Samples are transferred from the sput- 
tering chamber to the AFM in an environmental chamber filled with 
N2 in order to minimize the exposure to humidity and SKPM mea- 
surements are performed in dry and inert atmosphere (N2, relative 
humidity < 10%). 

Results 

The experiment consists of three steps. First of all, the back metal 
(Ag) is connected to the electrical ground and a SKPM scan is per- 
formed on top of metal surface (a portion of back metal surface is 
intentionally not covered with ZnO for this purpose). This is done to 
measure the metal work function. A calibrated tip (with known work 
function) is used for this purpose. Depending on the specific sample, 
small changes in value (±0.1 eV) have been measured. No depend- 
ency of Ag work function on relative humidity was found. After that, 
the same procedure is repeated with another tip in order to measure 
its work function. By doing so, we are able to calibrate each tip. 
Secondly, the tip is connected to the ground and the top surface of 
the ZnO sample is scanned in tapping mode. As previously men- 
tioned, this procedure is done in order to discharge the surface before 
starting the SKPM measurement. Finally, the same portion of the 
surface is scanned with SKPM. A positive ac voltage V^c is applied to 
the tip. In this case, electrons (majority carriers in n-doped semi- 
conductors) are attracted to the top surface of the sample. Thus, the 
semiconductor is in the accumulation regime and the metal/semi- 
conductor system can be treated in a first approximation as a metal/ 
metal system, allowing for the determination of the Vcpd^^- 



Before starting the experiments on ZnO, SKPM measurements are 
carried out on silicon and germanium samples with known charge 
carrier concentration. This is done to vahdate the model. Samples are 
etched with HF in order to remove their native oxide. Doping 
concentrations in accordance with nominal values are recovered, 
confirming the reliability of this technique. After that, the same 
experiments are performed on ZnO. 

An AFM topography image of the surface of a ZnO thin film with 
the corresponding contact potential Vcpd image, acquired at a dis- 
tance of 20 nm, is shown in Fig. 3. Fig. 3a shows the typical poly- 
crystalline nature of ZnO sputtered thin film, with small grains of 
nanometric dimension (tens of nanometers). The contact potential 
Vcpd image in Fig. 3b does not show the same morphological features 
reported in Fig. 3a. The Vcpd value appears almost uniform, varying 
within 10% of the mean value. Similar Vcpd values with a deviation in 
the order of 5% are recorded after scanning the same portion of the 
surface for more than 1 hour. This result confirms that no charging 
phenomena are taking place. The same experiment was repeated one 
week after the first measurement, keeping in the meanwhile the 
sample in inert atmosphere. No sensible variation in the Vcpd was 
measured. 

In order to relate the Vcpd with the effective doping concentration 
of thin film ZnO, all the parameters in equation (6) have to be 
recovered. In Table 1, the value of Nc and x for room temperature 
RF sputtered ZnO as well as the value of (j) for Ag are reported. While 
the first two are taken from^^, the value of ^ for Ag is measured for 
each sample as previously pointed out. 

As long as the sample temperature is kept constant, these para- 
meters do not change with thickness and deposition conditions'^. 
Fig. 4a shows a two-dimensional map of the effective doping con- 
centration in a 150 nm thick ZnO layer. In Fig. 4b the charge carrier 



Table 1 | ZnO and Ag parameters extracted from experimental 
results 

Quantity Value and Reference 

ZnO effective density of states in the conduction 3.7- 1 0^^ €171-^ 
band Nq 

ZnO electron affinity x '4.5 eV^^ 

Ag work function 0 4.7 ±0.1 V 
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Figure 4 | (a) Two-dimensional map of the charge carrier density in a 150 nm thick ZnO film, (b) doping concentration versus layer thickness for a set of 
ZnO thin films (from 10 to 300 nm). All the samples are deposited under the same conditions (deposition rate 0.1 A/sec, Ar pressure 0.8 mTorr). The 
effective doping concentration does not correlate with the layer thickness, except for the 10 nm layer. For each sample, three different regions have been 
scanned. The standard deviation of the value of the effective doping concentration is reported. 



density is plotted versus the layer thickness (from 10 to 300 nm). 
Each sample is deposited under the same conditions (deposition rate 
0.1 A/sec, Ar pressure 0.8 mTorr). The average value as well as the 
standard deviation of the contact potential difference Vcpd measured 
on all the samples is reported in the supplementary material. The 
uncertainty in the quantification of the free carrier concentration of 
this technique is given by the calibration process i.e. the quantifica- 
tion of the Vcpd between the tip and the Ag surface. 

Doping concentrations in the order of 10^^ cm"^ are recovered for all 
the samples. This value is in agreement with Four Probe measurements 
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performed on polycrystalline ZnO samples deposited under the same 
conditions on an insulating substrate (resistivity p = 0.7-1.8 Qcm, 
which, given a measured electron mobility of 31 ± 3 cm^/Vs, gives a 
doping concentration of ~1- 10^^-3- 10^^ cm"^). Furthermore, the 
doping concentration does not correlate with layer thickness. This 
is vaUd for layers with a thickness larger than the electron mean free 
path in ZnO (close to 20 nm)^°. Below this limit, the proposed model 
is not valid since quantum effects become relevant and the transport 
regime is not diffusive anymore^\ Therefore it cannot be applied to 
quantify the doping concentration in the material. 
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Figure 5 | Two-dimensional map of the charge carrier density in ZnO thin film (thickness 100 nm, deposition rate 0.2 A/sec, Ar pressure 10 mTorr). 
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In order to further validate the technique, a thin film ZnO with a 
thickness of 100 nm, deposited with a higher deposition rate (0.2 A/ 
sec) and with higher Ar pressure (10 mTorr), is measured. It has 
been experimentally demonstrated that an increase in the deposition 
rate as well as in the Ar pressure causes the doping carrier concen- 
tration in ZnO thin films to increase^^. In this case, we measure an 
increase in the work function difference (approximately 200 mV, 
from an average value of 445 mV to 647 mV) between the metal 
and the ZnO with respect to the previous case. This difference is 
attributed to a change in the charge carrier concentration because 
no considerable difference in mean surface roughness (from 4 to 
4.2 nm) as well as in electron mobility (from 31 to 40 cm^/Vs) was 
measured. Fig. 5 shows the sample charge concentration image. An 
average charge carrier density in the order of 10^^ cm"^ is calculated. 
Therefore, an increase of two orders of magnitude in the doping 
concentration is achieved by changing the deposition parameters. 
Being in accordance with experimental results reported in literature^^ 
and with our Four Probe measurements (p = 0.9-10"^ Q-cm), these 
results confirm the ability of the technique to measure the doping 
concentration of ZnO thin layers deposited under different 
conditions. 

Conclusions 

We probed the charge carrier density in ZnO thin films with high 
spatial resolution by using SKPM. A simple electronic energy model 
is used to correlate the work function of the sample with the effective 
doping concentration in the ZnO thin films. Limitations of this tech- 
nique (stray capacitances from tip geometry, water layer presence, 
surface charges and adsorbates) are discussed and possible solutions 
are proposed. The technique is experimentally tested on a set of 
samples deposited by RF sputtering under different conditions. 
Finally, we showed that results inferred with this technique are in 
accordance with carrier concentration expected for intrinsic ZnO 
thin films deposited with two different deposition rate and Ar pres- 
sure and obtained with Four Probe and Hall measurements. 

Methods 

SKPM measurements are performed on ZnO layers deposited by Radio Frequency 
(RF) sputtering on top of Ag. A set of samples with different thickness (from 10 to 
300 nm approximately) is analyzed. The samples are deposited and their surface is 
pretreated (Ar"^ bombardment at 250°C), annealed at 120°C in a vacuum chamber 
and immediately measured. SKPM measurements are performed in inert atmosphere 
(N2) with a Cypher (Asylum Research, USA). A Si cantilever coated with Ti-Ir 
(ASYELEC-01) is used. A cantilever resonant frequency coq of ~71 KHz, Q factor of 
130 and a spring constant A; of 1.5 N/m were measured using a thermal noise 
method^'. 

The average value as well as the standard deviation of the contact potential dif- 
ference VcpD measured on all thin film ZnO samples are reported in the supporting 
information. 
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